Effects of Strong Interactions in a Half Metallic Magnet: 
a Determinant Quantum Monte Carlo Study 
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Understanding the effects of electron-electron interactions in half metallic magnets (HMs), which 
have band structures with one gapped spin channel and one metallic channel, poses fundamental 
theoretical issues as well as having importance for their potential applications. Here we use 
determinant quantum Monte Carlo to study the impacts of an on-site Hubbard interaction U, 
finite temperature, and an external (Zeeman) magnetic field on a bilayer tight-binding model which 
is a half-metal in the absence of interactions, by calculating the spectral density, conductivity, spin 
polarization of carriers, and local magnetic properties. We quantify the effect of U on the degree 
of thermal depolarization, and follow relative band shifts and monitor when significant gap states 
appear, each of which can degrade the HM character. For this model, Zeeman coupling induces, 
at fixed particle number, two successive transitions: compensated half-metal with spin-down band 
gap — ► metallic ferromagnet — > saturated ferromagnetic insulator. However, over much of the more 
relevant parameter regime, the half-metallic properties are rather robust to U. 

PACS numbers: 71.10.Fd, 71.30.H-h. 02.70.Uu 



I. INTRODUCTION AND BACKGROUND 

As a family of promising candidates for spintronics 
application, half-metallic (HM) magnetic materials have 
attracted much interest over the past two decade*^, 
de Groot et al£- discovered HMs computationally and 
popularized the HM possibility especially in Heuslcr 
and half-Heusler compounds. The distinguishing 
characteristic of a HM is that at the mean field level the 
Fermi level for one spin direction falls within a gap for the 
other direction, resulting in 100% polarized conduction 
and obvious spintronics applications. Motivated by 
the unusual magneto-optical properties, the Heusler 
class of intcrmetallic materials provided the most 
promising realizations of HM character. Thermal 
fluctuations degrade spin alignment and thus destroy 
the ideal HM. In addition, real materials effects such 
as crystal imperfections can couple spins and degrade 
the polarization; spin-orbit coupling destroys the true 
HM. These phenomena can be minor at low temperatures 
compared with the Curie temperature (which ranges 500- 
730 K in the half Heusler family NiMnZ (Z=Co, Pd, 
Pt, Sb) i and can be higher), but because proposed 
applications do not rely on 100% polarization, HMs 
remain viable for near-future electronic devices. 

A variety of experimental techniques, including 
positron annihilation, spin-resolved photoemission 
spectroscopy, Andreev reflection, and nuclear magnetic 
resonance, have been employed to assess the character 
and polarization level of proposed HMs, often with 
less than definitive results. Comparision is made 
with density functional theory (DFT) based electronic 
structure calculations, which still play the dominant 
role in the specification of HMs and in the search for 
additional HMs and for the even more exotic half metallic 
antifcrromagnets^— . HMAFs, better characterized as 
compensated HMs, have zero macroscopic moment so 



they can provide additional functionalities. Known or 
strongly anticipated HMs now span a diverse collection 
of materials with different chemical and physical 
properties. 

An important question about HMs, not yet well 
clarified, is the impact of dynamic interactions (which 
lie beyond DFT methods) on the character and the 
survival of the single-spin-channcl gap that defines the 
HM phase. One of the most fundamental consequences of 
repulsive on-site interactions is local moment formation 
and dynamics, the study of whose effects in metals has 
a long history (Anderson impurities, Kondo systems). 
Local moments in insulators are also rather well studied 
and present a largely distinct set of issues (gap states 
versus band resonances, magnetic activity)^. HMs, 
especially in oxides, bring in all of these issues: although 
one spin is gapped, it is not electrically an insulator 
since there is metallic screening, but the two spin 
channels are fundamentally distinct. There have now 
been a number of studies, in particular by Katsnelson, 
Chioncel, and collaborators, that indicate degradation of 
the HM gap - sometimes severely - by electron-electron 
interactions. One general picture is that interactions, 
typically pictured in terms of a magnetic polaron formed 
by charge carrier-magnon binding, lead to tailing of 
spectral density into the gap and sometimes to mid-gap 
states that would substantially degrade performance as 
a HM device component. 

Another issue, not explicitly addressed to any great 
extent in treatments of interactions, is the distinction 
between the two types of HMs. The study of the fully 
polarized ferromagnet, also referred to as a saturated 
ferromagnet, in which the minority spin states are empty, 
extends back to the seminal work of Edwards and Hertz £ 
Cr02 is the simplest example of this: high spin S=l 
Cr 4+ has no occupied minority states^. The "gap" 
in Cr02 is several cV wide, between occupied O 2p 
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states and the unoccupied minority 3d bands. For 
low energy considerations, dc transport for example, 
minority spin states are out of the picture. The more 
common HM involves d states on either side of the 
gap as well as the continuum in the metallic channel; 
the minority charge excitations are typically narrowly 
gapped; in intermetallics a commonly ooccurring value 
is ^0.5 eV and even in oxides is rarely much greater. In 
the saturated case, there are simply no minority charge 
excitations to consider. 

There have been several studies based on model 
Hamiltoniana«i"ii. In particular, work based on 
the s-d exchange model predicted the existence of 
nonquasiparticle (NQP) density of states near the Fermi 
level, which (mathematically) arise from the branch 
cut of the self-energy arising from electron-magnon 
interactions^. Chioncel et.al. illustrated the coincidence 
between NQP states and the imaginary part of local 
spin-flip susceptibility in the framework of the single- 
band Hubbard Hamiltonian and DMFT approach (see 
below) where saturated ferromagnctism is stabilized by 
the additional magnetic spin splitting mimicking the 
local Hund ruleii. An effective spin Hamiltonian was 
derived to account for the temperature and disorder 
dependence of the magnetic properties of half-metallic 
double perovskitesi^. Remarkably, Kondo screening 
was recently shown to stabilize ferromagnetic order and 
further result in a half metallic phase with minority-spin 
gap in the Kondo lattice model with antiferromagnetic 
coupling^. 

The LDA+DMFT (local density approximation plus 
dynamical mean-field theory) approach that is becoming 
widely used to treat strong interactions is not based 
on a model Hamiltonian, instead using LDA results 
as the non-interacting system together with a self- 
interaction correction. This approach has been applied 
by Chioncel and coworker s) 11,14 ' 15 to evaluate the effect 
of interactions on HMs. One of the most well studied 
intermetallic HMs, NiMnSbjiS was the first application 
of this combined technique* 1 ^ The calculated spectral 
density contained non-quasiparticlc (NQP) states within 
the minority gap, but above, rather than pinned to, the 
Fermi level, allowing it to survive as a HM. However, 
further investigations of NiMnSb and other Heusler alloys 
show that the magnetic moment per formula unit, the 
NQP spectral weight, and the total DOS are insensitive 
or only weakly sensitive*! 1 , to correlation effects. In 
striking contrast, correlation effects were found to play 
a vital role in zincblcndc VAs, which is calculated 
to be a ferromagnetic semiconductor (i.e. gapped in 
both spin channels) within LSDA or the generalized 
gradient approximation GGA, but predicted to be a 
half-metal ferromagnet due to band shifts produced 
by LSDA+DMFT-il Strong correlation effects are also 
obtained in magnetite, which is uncommon in HMs 
in having a majority-spin gapJ£ In the full Heusler 
compound M^VAl which is HM within LSDA- 1 ^ an 
LSDA+DMFT treatment of local interaction led to NQP 



states within the gap but below the Fermi level-i 5 - which 
would not degrade spintronics-related properties. At 
present the impact of interactions appears to be highly 
material specific. However, the models and the treatment 
of the interactions (always approximate in some way) has 
varied widely, so few questions are truly settled. 

The DMFT approach, which treats on-site interactions 
and dynamics in detail, has the shortcoming of neglecting 
intersite correlations. Hence spinwaves, or even short 
range spin order, that still contain strong correlations 
between neighbors, arc replaced by identical but 
uncorrelated spin fluctuations on each interacting site. It 
is unclear to what degree the misrepresentation of these 
excitations may affect the character of the interacting 
spectral density. 

In this paper we investigate the survival (or not) and 
character of HM phases based on a bilayer Hubbard 
model with unequal interlayer hopping for the two spin 
species. This model allows a substantial parameter range 
in which the non-interacting density of states (DOS) has 
a gap in only one spin channel, and will be described 
in detail in the next section. We note here that the 
layer index can be regarded equivalcntly as an orbital 
or band index, so that the bilayer Hubbard Hamiltonian 
provides a useful pedagogical link between single orbital 
and multiorbital models. 

In the case when the spin up and spin down band 
structures are identical (i.e. in the absence of underlying 
magnetic order), tuning of the interlayer hybridization 
in such bilayer models has been demonstrated to drive a 
variety of quantum phase transitions^. For example, at 
half filling, the ground state can have antiferromagnetic 
long-range order for small interlayer (interband) hopping 
t_i_, and enter a disordered valence bond phase with 
singlet correlations between electrons on two layers, for 
large t±. Likewise, the system can evolve through Mott 
insulating transitions^— as t± is altered. In the doped 
system, there is a topological reconstruction of the Fermi 
surface, which modifies the spin fluctuations and changes 
the superconducting gap symmetry^. Adopting spin- 
asymmetric interlayer hopping, our model introduces 
new avenues of behavior to be illuminated. 



II. MODEL AND METHODOLOGY 

A. Previous work 

Previous studies of HMs in the single-band Hubbard 
model focusedii mostly on the limiting case of saturated 
ferromagnctism achieved through an external Zeeman 
magnetic field B. With the underlying up and down 
spin bands being degenerate, the effective spin chemical 
potentials /v = /x ± B are chosen to depopulate one of 
the species. An alternate way to achieve a half-metal, 
one in which the two spins species can still have the 
same population so that the polarization is zero, is to 
alter the band structure so that a gap opens for just one 
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of the species (which we will choose to be the "down" 
spins). For example, in a one band tight-binding model 
on a bipartite lattice with band e(k), one can incorporate 
an additional staggered potential Vj<j = (— 1) J V a , where 
the (— 1) J has opposite sign on the two sublattices. This 
alternating potential mixes momentum states k and 
k + 7r and opens up a gap in the dispersion relation 
E(k) = ±v/e(k) 2 + V; 2 . By choosing V f = and V± ^ 0, 
for an appropriate choice of Fermi level, the down spin 
species is insulating while the up species is metallic. 

Such a staggered potential, however, couples strongly 
to antiferromagnetism since it provides a one body energy 
which favors an oscillating down spin density on the 
two sublattices. The dominant nature of the resulting 
magnetic response might obscure the determination of 
the effect of U on the transport properties. For example, 
for a half-filled square lattice Hubbard model which has 
a divergent antiferromagnetic susceptibility xo( 7r 7 7r ) ) the 
additional staggered potential immediately produces a 
state with long range antiferromagnetic order (LRAFO) 
which is 'trivial' in the sense that it does not arise from 
a spontaneous breaking of symmetry, but rather from 
the externally imposed potential. This LRAFO opens a 
Slater gap in the initially metallic up spin spectrum, so 
that U immediately, but in some sense trivially, destroys 
HM behavior. 



B. Spin-asymmetric Hubbard model 

We avoid this confusing aspect by considering instead 
the slightly more complex case of a two layer (or, 
equivalently, two orbital) square lattice Hubbard model 
with spin-dependent inter-layer (inter-orbital) hopping, 

H = - t (4m* C jma + h - C ) ~ ~ aB ) n ™° 

(ij)ma ima 

ia 

+ U^2(n imt - -)(n imi --). 

im 

(i) 

Here the additional index m = 1,2 labels two layers 
(orbitals) while i,j are site indices and a is the spin. 
The first term is an intralayer nearest-neighbor hopping. 
We consider a square lattice with intralayer hopping 
t = 1 to set the energy scale. t± a is a spin-dependent 
interlayer (interorbital) hybridization, and U is an on- 
site repulsion. The terms coupling to the density are a 
spin-independent chemical potential fi and a Zeeman field 
B. The repulsive on-site interaction term is written in 
particle-hole symmetric (PHS) form so that at /i = B = 
the system is half-filled (for each spin species), even if 

The nonintcracting limit U — has two bands for each 



spin, 

e^(k) = — t± a — 2t(cos k x + cos k y ). 

ej(k) = +t± a — 2t(cosk x + cosk y ). (2) 

For t± a < At these two bands overlap, yielding 
metallic behavior. However, for t^ a > At, Eq. [5] 
characterizes a band insulator with gap 2(t± a — At). 
This Hamiltonian, and choice of U = band stucture, 
represents a half-metal without polarization, and avoids 
the externally imposed antiferromagnetism which would 
arise from a staggered potential. The price paid is 
the introduction of the extra layer (orbital) degree of 
freedom m. The Hubbard Hamiltonian on a diamond 
chain is a one dimensional analog in which the interplay 
of perpendicular hopping and U can create a correlation- 
induced half metal for certain fillings 25 . 



C. Underlying magnetic order 

The manner in which magnetic order (a fundamental 
necessity for a HM) is built into the model may affect 
behavior. In our model, magnetic order is implicit in 
the spin-dependence of t± a but is otherwise unspecified. 
A signature aspect of our model is that for U=0 it is 
a specific realization of the schematic symmetric HMAF 
DOS, presented for example in Fig. 1 of Ref. (2(|. In 
their DMFT study of a fully polarized Bethe lattice 
Hubbard model, Chioncel et.al. used a Zeeman field to 
mimic Hund's rule coupling, thereby splitting the two 
spin directions until the minority band was empty- 
In LSDA+DMFT studies of suspected HM materials 
(NiMnAs, FeMnSb, VAs), Chioncel and coworkersiiii^ 
based their dynamical corrections on the spin-split LSDA 
bands. 

The goal of this paper is to determine the effect of the 
electronic correlation term U in the model Hamiltonian 
Eq. [TJ Specifically, we compute the spin- resolved spectral 
densities as functions of t± a , U, and temperature T. We 
also study the interplay of U and interlayer hopping t± a 
on the antiferromagnetic correlations. For simplicity we 
will set = and vary t±± so that spin up fcrmions are 
metallic and the spin down fcrmions can be tuned from 
metal to band insulator at t±i = At. 

At the PHS point /.t = B = 0, both spin 
species are half-filled (regardless of the values of U, t, t± 
or temperature T). This immediately implies the 
polarization is identically zero, so that our model system 
realizes the exotic half-metallic antiferromagnetism 
(HMAF)-. Although there has been to date no clear 
confirmation of novel HMAF materials, a variety of 
candidates have been proposed, including La2VCuOg^ 
as a likely member of the double perovskite system with 
two magnetic ions,— semiconductors doped with dilute 
magnetic ions,— and monolayer superlattices CrS/FeS 
and VS/CoS^. An unusual semi-Dirac half-semimetal 
arises in untrathin VO2 film s 30 ' 31 . Our model's bands 
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closely resemble the simplest realization of a HMAF as 
arising from two exchange-split ions that are antialigned, 
shown schematically and discussed in Ref. Q. After 
considering the fundamental PHS case, we will also 
present some results for non-zero Zeeman field B, 
where the system can be expected to acquire nonzero 
polarization. 



D. DQMC procedures 

We treat the interaction using the determinant 
quantum Monte Carlo (DQMC) technique 3 ^. DQMC 
is a numerically exact approach to solve interacting tight 
binding electron Hamiltonians like Eq. [TJ In comparison 
with DMFT, DQMC has the advantage of being able 
easily to incorporate and measure magnetic, charge, 
and pairing correlations between different spatial sites. 
On the other hand, DQMC has the drawback of being 
formulated on finite spatial lattices so that finite size 
effects must be assessed. DQMC also is limited by the 
fcrmion "sign problem"— (much more so than single 
site DMFT) which prevents the acquisition of data at 
low temperatures. Most of the results presented in 
this paper will be for two 8x8 layers and inverse 
temperatures ft = 6. We will show that the sign problem 
is somewhat alleviated for Zeeman field B ^ so that 
lower temperatures can be reached. 



E. Properties to be studied 

In order to distinguish metal from insulator, and see 
the effect of U on the half-metallicity, we will examine the 
single-particle density of states, N a (ui). This quantity 
is obtained by an analytic continuation of the local 
imaginary-time dependent Green's function G a (r) = 

E( Tc ja( r ) C jv(°))' that is, by inverting, 
j 



G a ( T ) 



du- 



-j8w _|- 1 



(3) 



using the maximum entropy method 3 ^. Our focus will be 
on the density of states at the Fermi surface N a {w = 0) 
to see if HM behavior survives at non-zero U. 

Although the system has no net polarizaton at B = 
fi = from the viewpoint of total up and down 
occupations, the distinction between spin directions 
induces a polarization of the conduction electrons, which 
is quantified by, 



P{E F ) 



N t (E F ) - N^Ep) 
N t {E F ) + N^Ep) 



(4) 



This quantity has been the focus of much experimental 
work, since HMs (at T=0) have 100% polarization 
which is what makes them so attractive for spintronics 
applications. 



We also study the dc electrical conductivity Ode, 
which is extracted from the current-current correlation 
function, 



A x:E (k,T) 



(5) 



where ^(1,0) = iiE ff (4fi,A ~ c L c \+x,*)- The 
conductivity is obtained using the approximate form of 
the fluctuation-dissipation theore m 36 ' 37 , valid at large j3, 

A ra (k = 0,T = /3/2) = Trovfc/0 2 (6) 

The magnetic magnetic structure factor 

S W = Jj £ e4q ' (1 ~ j) (("'t - ^)(njt - n n )) (7) 



ij 



is of interest as well. In the ordered phase the spin 
correlations ((«if — iij.Xnjt — n il)) are non-zero for large 
1 j and 5(q) grows linearly with the lattice size N at 
the appropriate ordering wave vector q. 



III. CHOICE OF PARAMETER RANGES 

Most HM magnets investigated to date are at most 
moderately correlated intermetallic compounds^ and 
DFT treatments may be quite realistic. For instance, 
photocmission experiments and resonant x-ray scattering 
has led to the estimate of Hubbard interaction U = 2 
eV for NiMnSb^, while a Wannier orbital analysis has 
shown that the bandwidth of the NiMnSb bands crossing 
E F is W ~ 4 eV. The ratio U/W ~ 1/2 puts this HM 
compound in the weakly to moderately correlated regime. 
In our model system, the non-interacting bandwidth is 
W = 8t at t± =0, which suggests U/t = - 4 is the 
relevant range to study. 

Correlation effects of this size can fundamentally 
change the physics of related tight-binding Hamiltonians, 
e.g. opening a correlation gap or inducing AFLRO and 
a Slater gap; see for example Ref. [39| and references 
therein. The near-neighbor square lattice tight binding 
model is unstable to arbitrarily small U, although this 
sensitivity is a result of the van Hove singularity (vHs) 
and perfect nesting of that model. The majority channel 
in our model however retains these same features (at 
fx = only, of course); in the minority channel the 
vHs are shifted to cither side of ^t=0. Even with one 
remaining vHs, we will show below that interactions 
play a much less dramatic role when one spin channel 
is gapped. 

There is, after all, new physics in a HM as well as the 
new phenomena discussed in Sec. I: the new energy scale 
that is given by the spin-down gap A = 2(tj_ 4. — At). 
Down-spin charge excitations are gapped by A, which 
make interactions further differentiate the spins, for 
example by enabling sharp 'impurity' states (only) within 
the spin-down gap. Spin-flip interactions are gapped by 
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FIG. 1: Half metallic behavior persists in the presence of 
interactions: the down species density of states is virtually 
unaltered at U/t = 4 from its U/t = value. Here t±^ = and 
t±i — 5 so that the down species (only) has a noninteracting 
band gap A = 2f. The up species is metallic with a density 
of states at the Fermi surface that is relatively insensitive to 
U. The lattice size N = 8 x 8 in each layer and the inverse 
temperature (8 = 5. The chemical potential and Zeeman fields 
/i = B — 0. However, due to the particle-hole symmetry at 
half filling, there is no signature associated with the NQP 
states. 

the separation between [i (the up-spin Fermi level) and 
the nearest down-spin band edge, which is A/2 at ^=0 in 
our model. Kondo processes vanish due to the absence of 
low energy spin-flip excitations. Perhaps more relevantly, 
the vHs- and nesting-driven spin density wave instability 
obviously will be quenched when there is a minority gap. 
We show it is also suppressed in the ferromagnetic metal 
phase where the minority bands are split but not gapped. 

IV. EFFECT OF ON-SITE INTERACTION 

The moderating of strong interation effects is seen 
in the density of states, shown in Fig. [TJ where 
the insulating species N(u>) at U/t = 4 is virtually 
indistinguishable from its U — form. Here we have 
set t±f = and t±± = 5t and [i = B = so that 
the down spin species has a noninteracting band gap 
A = 2(t±i — At) = 2t and the up species is metallic. 
The gap magnitude and tailing of states into the gap is 
unaltered up to U=4. The up (metallic) species Nf(Ep) 
is also constant from U/t = to U/t = 4 to within the 
accuracy of the maximum entropy inversion of Eq.[3J The 
peaks in Nf(ui) away from ui = arc affected by the finite 
size effect of the 8x8 lattice. 

One distinguishing characteristic of interactions in 
HMs extensively discussed in previous studies is the 
appearance of the non-quasiparticle (NQP) states within 
the spectral gap - the magnetic polaron effect. These 
NQP states in previous studies have (1) appeared, at u = 
0, above w = 0, and below uj = 0, and (2) also possibly 
not appeared. In the generic picture, the density of NQP 
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FIG. 2: The dc conductivity for /3 — 6 is shown as a function 
of U. For the half metal t±± = 5t, U has very little effect on 
CTdc, which is consistent with the invariance of the density of 
states in Fig. [1] In a fully metallic phase t±i — there is a 
clear decrease of <Tdc with U due to the additional electron- 
electron scattering which can occur when both species have 
nonzero density of states at the Fermi surface. 
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FIG. 3: Staggered static structure factor Saf through the 
metal-HM transition. Increasing interlayer coupling reduces 
the increase with increasing U until, for half metallic phase 
t±i = 5, U has little effect on the (small) magnetic 
coupling. In symmetric bilayers, large interlayer promote 
singlet formation, which are suppressed here by the HM gap. 

states vanishes at the Fermi level (w = 0) but increases 
toward an energy scale of the order of the magnon 
frequency, leading to an asymmetry of spectral function^ 
Assisted perhaps by the particle-hole symmetry in our 
model (implying the symmetry of spectral function) at 
half filling, NQP states may be inhibited from appearing, 
and there is no signature of NQP states in Fig. [TJ The 
implicit nature of magnetic order in our model may also 
play a role, but other treatments have also incorporated 
some implicit origin of the magnetic order in a HM. 

The HM character can be expected to become more or 
less evident in a property specific manner, and we now 
describe a few examples. The relatively minor effect of 
on-site interaction U is further evidenced by the behavior 
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of the dc- conductivity in Fig. where the interlayer 
coupling is varied to move the system through the 
metallic phase tj_j_=2, 3, through the transition t_n.=4, 
to the HM phase t±^=5. At t±^=2 the conductivity is U- 
dependent because N(lu — > 0) in the metallic phase has 
the underlying van Hove singularity there. As U increases 
the dependence weakens, and at t±i=4 and 5, which 
is crossing over into the HM phase, Ode decreases and 
becomes independent of the interaction strength: in the 
HM phase t±i = 5, <7dc falls by less than ten percent as 
U increases from U = 1 to U = 3, and must be controlled 
solely by spin-up processes. 

Fig. [3] demonstrates the correlation effects on the 
staggered magnetic static structure factor, again for the 
progression from metal to HM t±±=2-5. The increase 
with U at small interlayer coupling decreases as the HM 
phase is approached and entered, and the variation of 
5af with U in the HM phase t±i — 5 is quite small. It 
is unlikely that antifcrromagnctic long-range order will 
arise at lower temperatures, and it is known that, in 
the case of spin-independent interlayer hybridization, t± 
drives a competing singlet ground state, with a quantum 
critical point at tj_ ~ 1.6 above which the ground state 
of the bilayer model no longer has AFLRO in its ground 
state. 

In summary, our simple bilayer model of an 
unpolarizcd HM indicates the robustness of the HM 
phase to interactions. This is consistent with some 
previous work, which concludes that the half-metallic 
properties of several materials such as NiMnSb and other 
Heusler alloys are insensitive or only weakly sensitive 
to correlation effect*^. It is not apparent whether the 
explicit particle-hole symmetry at half filling and implicit 
nature of magnetic order in our model has any substantial 
effect, but whatever the reason there is no signature here 
associated with the states in the gap. 



V. EFFECT OF TEMPERATURE 

The temperature dependence of the half-metallic 
properties and their stability against finite-temperature 
spin excitations are crucial for practical applications. As 
mentioned before, one crucial effect is the depolarization 
caused by finite-temperature. With the separate 
spin bands fixed in our model, the effect of reduced 
magnetization upon approaching the Curie temperature 
is not included, but we can study temperature effects "at 
fixed magnetization." 

Fig. @] illustrates the temperature effects on the density 
of states in the half-metallic phase. While substantial 
rearrangements of spectral weight occur for spin up, the 
primary effect for spin down is an increased tailing of 
states across the band edges and into the gap as T 
increases. Only when /3 < 3 is the half-metallic feature 
destroyed, i.e. N\(uj — > 0) becomes appreciably non- 
zero. This corresponds to a temperature T ~ t/3 ~ 
A/6 since the down spin band gap A = 2t. If we 
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FIG. 4: Main panel: Temperature effects on the density of 
states of half- metallic phase (U = 3). The inset shows N±(u = 
0) versus T. The half- metal gap in Ni(uj) appears to be robust 
for > 3 (T < t/3). 




T 



FIG. 5: The temperature dependence of the conduction 
electron spin polarization, where the non-interacting gap is 
equal to A = 2t (t=l). P(Ef) begins to turn downwards at 
T ~ A/6, with minor dependence on the value of U. 



phenomenologically introduce A(T) which vanishes at 
the Curie temperature, and suppose that the magnitude 
of the gap is the primary energy scale for this purpose, 
then we can infer destruction of HM character by 
interactions (and thermal broadening) around T* = 
A(T*)/6. 

We can also examine thermal effects by evaluating 
the conduction electron spin polarization P(Ep). The 
combined effects of finite-temperature and Hubbard 
interaction U are given in Fig. [SJ The polarization begins 
to deviate downwards from unity at T ~ t/3. The 
role of U on P(Ep) is negligible in the low temperature 
regime, but as T increases to intermediate temperatures 
U has a distinct depolarization effect. Thinking of our 
parameters as very roughly relevant to NiMnSb, t/3 is 
a very high temperature, one at which the underlying 
magnetic order (fixed in our model) physically has 
diminished drastically or vanished. Even well above 
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FIG. 6: Temperature dependence of conductivity, for U=2 
and 4. In the HM phase (t±± = 5) the conductivity (see Eq. 
6) is contributed only by the up spin carriers, which as noted 
earlier is not affected strongly by these values of U. 



room temperature thermal effects as well as effects due 
to interaction U up to 4 on P(Ep) are negligible. 

A separate diagnostic of the nature of half-metallic 
phase and its robustness against the Hubbard interaction 
U is provided by the conductivity, Fig. H3 evaluated from 
Eq. 6. Recall that in a half metal for these values 
of parameters, the conductivity arises only from the 
majority (ungapped) channel, and spin-flip scattering is 
frozen out at these temperatures leaving scattering only 
in the charge channel. This data supports the previous 
indications of the relatively small effects of interactions in 
our bilayer model of a half metal. The primary difference 
for U = 2 is the significantly larger conductivity when 
T falls below 0.2t, compared to the U = 4 trajectory. 
On the other hand, adc at higher temperatures hardly 
depends on U. 



VI. EFFECT OF ZEEMAN MAGNETIC FIELD 

All of results above are at fj, = B = which, by 
particle-hole symmetry, guarantees half filling of both 
spin species: (m m a) = 0.5. A Zeeman field B introduces 
a spin bias, but one of the signatures of a HM is that it 
is impervious to magnetic fields that are not too large, 
that is, as long as the gap persists and fj, does not cross 
a gap edge. This vanishing of the spin susceptibility 
is self-evident at the mean field level, where B merely 
shifts the relative positions of the up and down bands. 
If the chemical potential (finally determined by state 
filling in the metallic channel) does not cross either gap 
edge, there is no rcoccupation of states of either spin. 
In fact there is no change (up to a constant) in the real- 
space potential for either spin, so the (many-body) states 
themselves do not change. The energy changes only due 
to the trivial magnetic energy term —MB (M is the 
net magnetization, which is unchanging until fi crosses 
a band edge) which is zero in our HMAF-type model. 
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FIG. 7: Effects of Zeeman field B on the density of states, at 
fixed fj,. Here U = 3, /3 = 8, t±i = 5. The Zeeman field B 
shifts the spectral function of both spin species. Top panel: 
At weak fields (viz. less than half the gap), N(u) more or 
less rigidly shifts with B for both spin species, as in mean 
field. In the case of the insulating down spin electrons this 
shift eliminates insulating behavior as B approaches half the 
band gap 2(t± — At) and the magnetization begins to change. 
Bottom panel: At large magnetic field the up spin density of 
states is driven completely below the Fermi level /i + B so 
that now the up species is insulating. Meanwhile the down 
spin density is metallic because, with its larger bandwidth, 
the down spin Fermi level [i — B is not yet completely below 
the down bands. In this way one has a HM in which the 
metal/insulator roles of the two spin species is reversed. 



The interest here is in interacting systems, and this 
anticipated behavior of HMs has been supported by 
some rigorous results based on ground state many-body 
wavefunctions and spin density functional theor y 40 ' 41 . 
The theory provides for ranges of applied Zeeman fields 
for which the ground state, and therefore the spin-density 
matrix is unchanging, just as at the mean field level. 
The values of the applied fields, positive and negative, 
at which the magnetism changes thereby provides the 
gap edges, and their difference provides the (interacting) 
fixed particle number gap. 



In our calculations when t 



-it 



t±l and B 







there is no sign problem at any temperature. This is 
a consequence of a total correlation between the signs 



8 



of the up and down spin determinants, so that their 
product, the probability of the configuration, is always 
positive. Allowing p or B to become nonzero induces 
a sign problem so that, normally, simulations are much 
more challenging. Here, however, with tx_-\ ^ t±i the 
correlation between spin up and spin down determinants 
has already been broken, which is why our simulations in 
the earlier sections do not extend beyond /3 w 6. We can 
access a similar range here when B ^ 0. 

In mean-field level, the effect of the external Zeeman 
field is only to shift the spectral functions of both spin 
species rigidly in opposite directions. It is therefore 
natural to expect a transition from half-metallic phase 
to normal phase at \B C \ equal to half the band gap. 
We show that this transition survives in the presence 
of Hubbard interaction U, but with spectral weight 
redistribution which is not captured in MF in addition 
to renormalization of B c . 

Fig. [7J shows the effects of Zeeman field on the spectral 
density. In the top panel the weak magnetic field-induced 
shifting of the spectral weights of both spin species in 
the opposite directions is clear, though the shifts begin 
to deviate from being rigid. HM character survives 
to B c w0.4, compared to £> c =0.5 without interactions. 
This "magnetic field gap" is renormalized by ~ 20% at 
U=3. In the bottom panel, we demonstrate that a strong 
enough magnetic field (of order the half the bandwidth, 
B W/2) can induce a situation in which the up spin 
electron bands become completely filled (heuristically, 
the effective chemical potential p c ^ = p + B lies above 
the band), while the down spin electrons remain metallic: 
Mcff^ = —p+B still cuts across a region of nonzero N^(lo). 
Since the filled majority channel holds one electron, the 
occupation of some minority states reflects the fact that 
the total density has increased as the field is increased at 
fixed chemical potential p=0. At fixed particle density, 
a large field will indeed produce a filled majority band, 
in which the minority band must be empty (the non-zero 
spectral density must lie only above p—0). This state 
is a saturated ferromagnetic insulator. In this way, as 
B increases, the sequence of phases at fixed p proceed 
from HM with majority carriers — > metallic ferromagnet 
— > HM with minority carriers. This sequene converts, at 
fixed density, to HM with majority carriers — > metallic 
ferromagnet — > saturated ferromagnet insulator. 

Another reflection of the magnetic behavior, again 
at fixed p=0, is illustrated in Fig. [5J The external 
Zeeman magnetic field leads to an increasingly polarized 
lattice [top panel, P = (pf - pi)/(p\ + pi)), as p T grows 
monotonically at the expense of p^. At B ~ 4 the up 
bands are completely full while the down bands, which 
due to their nonzero have a larger bandwidth, are still 
crossed by the effective chemical potential. This appears 
to be re-entry into a (polarized) HM phase. However, the 
spectra in Fig. [7] and the data of Fig. [5] were obtained at 
constant p, and the B-induced increase in density has bc 
recognized but docs not affect this picture. 




0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4 
B 

FIG. 8: Up, down, and total density as function of Zeeman 
field at fixed /i = and increasing B. The system starts at 
half-filling p^ = pi = 0.5 at B = and becomes increasingly 
spin polarized as B increases. Ultimately at B ~ 4 equal 
to half the band- width, the up bands are completely filled. 
This value of B is not quite sufficient to empty the down 
band owing to its slightly larger band-width. The polarization 
is P = (p-\ — Pt)/(Pt + Pi)' The parameters are chosen as 
U = 3,@ = 8. 



VII. SUMMARY 

We have used the numerically exact finite-temperature 
determinant Quantum Monte Carlo (DQMC) method to 
study the effect of strong interaction induced correlations 
on half-metallic behavior in a multiband Hubbard 
Hamiltonian. Our model consists of a bilaycr square 
lattice, or, alternately viewed, a two orbital system, 
with spin-dependent interlaycr hybridization chosen to 
induce a band gap in only one spin species. This is 
a model appropriate to half metallic antifcrromagncts, 
since the lattice has the same number of spin up and 
spin down electrons, but only the latter have a gapped 
non-interacting spectrum. 

By investigating the influence of an on-site Hubbard 
interaction U, finite temperature T and external Zeeman 
magnetic field B, we find that the half-metallic properties 
arc not particularly sensitive to the interaction U, up 
to values equal to about half the bandwidth which 
seem appropriate to intermctallic half metals. Finite- 
temperature effects depolarize the conduction electron 
states only at T > i/3 of the gap, with a degree of 
depolarization which depends weakly on U. A very 
large Zeeman magnetic field drives the system (at fixed 
particle number) from half metal to metallic ferromagnet, 
and finally to a ferromagnetic insulating phase when 
the minority spectral density is Zeeman split completely 
above the majority spectrum. 

An interesting issue in half-metallic ferromagnetic 
materials that remains undecided is when non- 
quasiparticlc (NQP) states arising from electron-magnon 
interaction arise, and whether they are above, below, or 
spanning the chemical potential. The signature of such 
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states is the appearance of a resonance in the gapped 
channel. In contrast to several earlier studies based on 
different models and using a different treatment of the 
interaction, we have not seen any distinct characteristic 
features associated with NQP states in the DQMC 
studies of our model. The difference may be related 
to the fact that our bilayer system without the external 
magnetic field satisfies particle-hole symmetry so that Ef 
is always pinned to the center of the gap. The other 
possibility is that the sign problem has prevented us 
from reaching low enough temperature to observe the 
development of the resonance, or that the spin waves 
that may be required for a proper description of electron- 



magnon interactions arc not yet fully formed. It is 
possible that the effects of interactions are simply model 
and material dependent rather than universal. 
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